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Introduction. A method for solving the problem of an infinite 
plate on an elastic foundation is proposed. The plate is affected 
by a periodic load in the form of a force following an arbitrary 
closed path. The work objective is to develop a numerical meth- 
od for solving problems of the elasticity theory for bodies under 
a moving load. 

Materials and Methods. Given the periodicity of the load under 
consideration, it is decomposed in a Fourier series in a time in- 
terval whose length is equal to the load period. The solution to 
the original problem is constructed by superposition of the solu- 
tions to the problems corresponding to the load specified by the 
terms of the Fourier series described above. The final solution to 
the problem is presented as a segment of a series. In this case, 
each term corresponds to the solution of the problem of the im- 
pact on an infinite plate of a load distributed along a closed 
curve (the trajectory of the force motion). To find these solu- 
tions, the fundamental solution to the equation of vibration of an 
infinite plate lying on an elastic base is used. 

Research Results. A new method is proposed for solving prob- 
lems on the elasticity theory for bodies with a load following a 
closed path of arbitrary shape. The problem of an infinite plane 
along which a concentrated force moves at a constant speed is 
solved. It is determined that the trajectory of motion is a smooth 
closed curve consisting of circular arcs. The behavior of dis- 
placements and stresses near a moving force is considered. The 
energy propagation of the elastic waves is studied. For this pur- 
pose, the coordinates of the Umov — Poynting vector are calcu- 
lated. The effect of the force motion speed on the length of the 
Umov — Poynting vector is investigated. 

Discussion and Conclusions. The method is applicable when 
considering more complex objects (plates of complex shape, 
layered plates, viscoelastic plates). Its advantage 1s profitability 


since the known problem solutions are used to build the solution. 
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Beedenue. IIpennaraetca MeToy pellieHua 3a7a4u 0 OecKOHeuHON 
TulacTHHe, JIe*Kallled Ha yupyroM ocHoBaHUuH. Ha miacTuHy jen- 
CTByeT MepHoyMueckad Harpy3Ka B BUJe CHJIbI, WepeMelllaroleiica 
110 IIpOW3BOJIbHOH 3aMKHYTOM TpaekTopnn. Lew uccneqoBaHua 
— pa3pa0oTKa 4HCJICHHOrO MeTOJa pellleHHaA 3ajad Teopuu yripy- 
TOCTH JIA Tes, HaxOJAMIMXcA MOA evCTBHeM TIOJBYDKHOM 
Harpy3Ku. 

Mamepuaiei u memoool. YuviTbiBad TepHOAMYHOCTb paccMaTpu- 
BaeMOM Halpy3KH, OHA packJlaj[bIBaeTca B pa, Dypbe Ha BpeMeH- 
HOM OTpe3Ke, JJIMHa KOTOpOrO paBHa Mepwory Harpy3ku. Perrte- 
Hue UCXOAHONM 3aa4yH CTpoOwMTcA MOcpexCTBOM cynepio3snuuu 
pelleHuit 3aja4, COOTBETCTBYIOINMX Harpy3Ke, 3a/aBaeMon CJlara- 
€MBbIMM OMMCaHHOro BbilIe pxya Dypbe. OKOHYATEIIBHOE pellleHHe 
3aa4u MpeyIcTaBIAeTCA B BUe OTpe3Ka pala. Kaxkyoe ciaraemoe 
IIpH 39TOM COOTBETCTBYeT PelIICHHFO 3aja4H O BO3ZeMCTBHU Ha 
OeCKOHeYHY!O IvIacTHHy Harpy3KH, paciipeyesyIeHHOW MO 3aMKHy- 
TOM KpWBOM (TpaeKTopHu BYwKeHHA CHIBI). JIA HaxoxKeHHA 
9THX pellieHH NcHONB3yeTca (PyHTAaMeHTasIbHOe pellleHve ypaB- 
HeHuA KOeOaHHA OeCKOHeYHON MaCTHHBI, JieKallle Ha yipyroM 
OCHOBaHHH. 

Pezyiomamoi ucciedoeanua. [TpenyioxKeH HOBbIM MeTO pellleHna 
3alad TEOpuN yNpyrocTH JIA Tel C Harpy3KOH, JBYDKyIelica 10 
3AaMKHYTOM TpacKTOPHU MPOu3BOJIbHON PopMbI. PeiieHa 3ayaua oO 
Oe€CKOHEYHOH MWIOCKOCTU, TO KOTOPOHM C MOCTOAHHOM CKOPOCTbIO 
J{BMMKETCA COCpeqOTOUeHHad cua. OrpeyeseHo, 4TO TpaeKTOpuA 
TBWKEHUA TpescTaBsiaeT COOOM rylayKy}O 3AMKHYTYIO KPHBYIO, 
COCTOALLY1O H3 Yr OKpy2xKHOcTeH. PaccMoTpeH xapakTep H3MeHe- 
HUA TMepeMellieHHH MW HallpsKeHuM BOM3H JBYDKYyUelica CHIBI. 
V3y4eHo pacipoctpaHeHve 39HeEpruu yupyrux Bos. C 3TOM Ie- 
JIb}O BbINOJIHEHO BBIMMCJIeHHe KOOpyMHaT BeKTOpa YMoBa — 
Tlovtunra. UccieqoBaHo BIIMAHVWe CKOPOCTH JIBYDKeHHA CHJIbI Ha 
JyIMHy BeKTOpa YMosBa — IlovTuura. 

O6cyorcoenue u 3aKo4eHuaA. Mero, 1pHMeHHM HU IIpu paccMor- 
peHuU Oosiee CIOXKHBIX OOBECKTOB (TMIMTHI CO2KHOM (DOPMBI, CIIOU- 
CTbIe¢ IWIMTbI, BA3KOYIpyrnve WMTbI). Ero npenMyljecTBO — 9KO- 
HOMHUHOCTE, Tak Kak JIA MOCTPOeCHHA pellieHHA MCTIOJIb3y!OTCA 
yoKe M3BECTHBIe pellleHua 3ayjay. OKOHYAaTeIbHOe pellleHve BbIpa- 
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The final decision is expressed in a convenient form — as the 
sum of curvilinear integrals. The results obtained can be used in 
the road design process. Studying the energy propagation of 
elastic waves from moving vehicles will enable to evaluate the 
impact of these waves on buildings near the road. The wear of 
the pavement is estimated considering data on the behavior of 


displacements and stresses. 
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%KaeTCA B YIOOHOM BUyIe — KaK CyYMMa KPHBOJIMHeCHMHBIX HHTerpa- 
nos. [losyueHHble pe3yiIbTaTbI MOryT ObITb UCIOJIb30BaHBI B TIpo- 
yecce mpoekrupopaHua yopor. M3yyenve pacnmpocrpaHenua 
9HeEprHu YIpPyIx BOJIH OT JBMKYIIMXCA TPaHCIOpTHBIX CpeCTB 
TIOSBOJIMT OLCHHTb BO3eMCTBHe yKa3aHHbIxX BOJIH Ha CTPOeHHA, 
pacnosioxKeHHble BOIM3H WOporu. C y4eToOM J{aHHBIX O XapakTepe 
VW3MeHeHHA TepeMelleHHH UW HalipsxKeHHM OeCHMBaeTCA W3HOC 
JOPOXKHOLO NOKPBITHA. 
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Introduction. The study on dynamic phenomena caused by the action of a moving load is an urgent task that 
has application significance (for example, when solving transport development issues). In the papers devoted to this 
problem, various tasks with a moving load were considered. In particular, it is shown how a load moving in an infinite 
straight line at a constant speed acts on a half-plane or half-space (elastic isotropic, transversal isotropic, viscoelastic). 
In this case, when solving the problem, a moving coordinate system associated with a moving force is introduced; this 
enables to exclude time from the number of independent variables [1—5]. Some papers consider the action on an infinite 
plate or strip (elastic or viscoelastic) moving uniformly along a rectilinear load path. In this case, the same method of 
eliminating a temporary variable is used, or a quasistatic formulation of the problem is considered [6—12]. The tasks in 
which the length of the load path is finite and the path itself is a curved line are of main interest. In this case, finite ele- 
ment modeling of a moving load is frequently used [11—13]. In a number of papers, when solving such problems, varia- 
tional methods are used (in particular, the Rayleigh — Ritz method) [14-16] or a variation of the Galerkin method, 
which makes it possible to reduce the problem to ordinary differential equations. In this case, various objects of the ap- 
plication of a moving load are considered (plates, layered plates, viscoelastic plates, half-spaces — both isotropic and 
anisotropic) [17—19]. This paper presents a method that develops the ideas described in [20—23]. 


Materials and Methods. Consider an infinite plate lying on an elastic Winkler base, which is under the action 
of a normally applied force moving along a closed path. 
The problem is reduced to the integration of the equation of motion of a plate lying on an elastic Winkler base [14]: 


: P 
AW+0°0W+kW =~, (1) 
D 
where W is the plate deflection; D = Dev): E is Young's modulus; v is Poisson's ratio; h is the plate thickness; 
—V 


: h , KC. ; ; 
c= ; p is the density of the material; k = = ; K is the compliance coefficient of the elastic base; P is the concentrated 
force moving along a closed curve y with the constant speed a. 

Let us introduce the coordinate s counted from some fixed point of the curve y. Then the force P moving along the 
curve y with the velocity a will be described by the relation P = P(s —at). The function P(s — at) 1s periodic in t, with peri- 


L 
od 7T = — 
a 


, where L is the length of the curvey. 
Solution. Consider the steady state. We expand the function P(s—at) in the Fourier series in the variable ¢ on 


LoL ; 
the segment S05 | . In this case, the expansion coefficients appear as: 
a 2a 
Z 
72a Ziaknt/ 
Cc, = | P(s —at)e Ldt . 


Yo 


After changing the integration variable in the integral s —at = z , we obtain: 
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L ; : 
Ly -2iks/ a aes 
Cc. = | P(z)e Ldz =d_, ; 


a 
“7” 


; ; LL 
where d, are Fourier series expansion coefficients on a function segment 55 | of the function P(z). 








Then the moving load can be presented as a Fourier series: 


1 = 2ikn(s—at) 
P(s-at)=— ) dye Y, ; 
aa 


=—0O 


Given the linearity of the problem, its solution can be presented as: 
W = » d_ W, 5 (2) 
k=—0o 
where W, are plate deflections caused by the action of a vertical load whose distribution along the curve y is described by the 


2, ae pci 
functione 7 that varies in time according to the law e 


To determine W, , Wwe use the fundamental solution to the equation (1), which corresponds to 


P=85(x-X))S(V— oe" , THe o, = “en 


Using the limiting absorption principle and traditional methods for constructing solutions to differential equations, 
2 


‘ i ; : 60) 
we can obtain a fundamental solution to the equation (1), which at k > ss , has the form: 
C 





Wy (X, V3 X9> Vo) = papi (a,R)—Ky (aR) | ? 


j 2 a _jT 
where R= [(x-x y +(y-¥ y , oy le /, , = cel , O, =Ee i , Ky (z) is the Macdonald function. 
% 


2 
w oe oe oe 
Atk< — , the solution is given by: 
C 





 11$)(1R)~ Ko (18) |, 


L 


Ary” D 


2 
where y = 4 —k , HS” (yR) is the Hankel function. 
C 


7 2ikn(s—at) 
Then W, = po, (X, V, Xo (5), Vo(s Ne Y, ds . Using the well-known formulas of the thin plate theory and 
: 
the formula obtained from the above relations determining the deflection W (2), we can calculate the displacements w,, 


u, and the stresses 6,, 6, and 6,, at any point on the plate. 
For large k, it is necessary to calculate the integral of fast oscillating functions. For this, a quadrature formula 


based on replacing the weakly oscillating part of the integrand by a cubic spline was used, and the highly oscillating 


a 
factor e ’ was considered as a weight function [15]. 
Research Results. So, an infinite plate lies on an elastic Winkler base. A normal force acts on it moving along 


the trajectory shown in Fig. 1. 
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Fig. 1. Load motion trajectory 


The computation was performed at the following initial data: h = 0.25 m; s = 221 m/s; E = 232469°10 N/m2; v 
= 0.36. Fig. 2 shows the calculation results corresponding to K = 1.864 m “ and a = 125 m/s. 
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Fig. 2. Variation of displacements and stresses 


These graphs describe changes in displacements and stresses in the coordinate system associated with the point 
of application of the moving force P. The ¢ axis is directed tangentially to the motion path? , and the n axis is directed 


along the external normal to the region bounded by the path 7 (see Fig. 1) at z= VA , 
In this case, the displacement vector and stress tensor, respectively, were represented as: 
U =Ut-t+Un-n+Ww-k , S = St-7f +Sn-7m + Sin-(H +70), 
where & is normal to the plate. 


Fig. 2 shows the change along the ¢ axis of the displacement vector components Wt,Utt,Unt , the stress ten- 


sor Sit, Snt, Sint , and the variation of these values along the n axis— Wn,Utn,Unn and Stn, Snn, Stnn . 
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According to the calculations, the behavior of these values at all points of the trajectory remains unchanged. 
The computation also showed that with a change in the velocity of the load a in the range from 0 to 125 m/s, the com- 
ponents of the displacement vector and the stress tensor increased moderately (by 3-4%). 

To study the elastic energy propagation, the components of the Umov-Poynting vector e, =—S;,,-0,u; (Sj are 


the components of the stress tensor, u; are the coordinates of the displacement vector) whose direction indicates the 


direction of energy propagation, and the length describes the amount of energy transferred through a surface unit per- 
pendicular to this vector line per unit time. 
Fig. 3 shows the elastic energy propagation near a moving concentrated force (the position of the force is 


marked with a red asterisk. 





Fig. 3 Elastic energy propagation 


Discussion and Conclusions. The analysis of the results obtained shows the following: at the speed variation 
limits indicated above, the length of the Umov-Poynting vector is almost proportional to the load speed and a sound 
energy pattern near the force changes slightly during the movement. The behavior of the displacements and stresses 
calculated above at all points of the trajectory remains unchanged, and their values weakly depend on the load speed a 
when this velocity varies from 0 to 125 m/s. 

Application of the proposed method to an infinite plate lying on an elastic base does not exhaust its possibili- 
ties. It can be used when considering more complex objects (plates of complex shape, layered plates, viscoelastic 
plates). The considered method differs from the mentioned above in greater efficiency since it uses the known problem 
solving to construct the solution. The final decision is expressed in a convenient form — as a sum of curvilinear inte- 
grals. 

The results can be used in the road design process. Studying the elastic energy propagation from moving vehi- 
cles will provide evaluating the impact of these waves on buildings near the road. The wear of the pavement is estimat- 
ed considering data on the behavior of displacements and stresses 
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